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Abstract. 

We present the result of 28 hours of observations of the nearby starburst galaxy NGC 253 with the H.E.S.S. detector in 2003. 
We find no evidence for very high energy y-ray emission from this object. Gamma-ray emission above 400 GeV from NGC 253 
had been reported by the CANGAROO collaboration in 2002. From the H.E.S.S. data we derive upper limits on the flux above 
300 GeV of 1.9 x 10"'^ photons cm"^ s"' for a point-like source and 6.3 x 10"'" photons cm"^ s"' for a source of radius 0.5° 
as reported by CANGAROO, both at a confidence level of 99%. These upper limits are inconsistent with the spectrum reported 
by CANGAROO. The expected very high energy y-ray emission from this object is discussed in the framework of a galactic 
wind propagation model. 
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1. Introduction 

Starburst galaxies, such as the closest objects of this type, 
NGC 253 in the Southern and M 82 in the Northern 
Hemisphere, are known to exhibit enhanced and strongly lo- 
calised supernova explosion rates in so called 'starburst re- 
gions'. In such regions gas and photon densities are very high 
compared to the situation in normal star forming galaxies like 
the Milky Way. Given the paradigm of cosmic -ray (CR) ac- 
celeration in supernova remnants (SNR), substantially greater 
CR energy densities > 100 x Ef^' (where E^^^ ~ 1 eV cm"^ 
is the value in our galaxy) are expected in such regions, with 
corresponding fluxes for 7r°-decay y-rays alone in the range 
of lO-*^ < Fy(> 300 Me V) < 10 "^ pho tons cm^^ s"' for an 
object such as M 82 (Volk et al. [1989"). Later estimates for 
M82 (Akyuz et al. ll99U Volk et al.ll996l are approximately 
encompassed by this range, which translates into a y-ray flux of 
1.5 X 10-'^ < Fy(> lTeV)x(£^/lTeV)i i < 1.5 x lO^'^ pho- 
tons cm"^ s"' at energies Ey » 300 MeV assuming a differen- 
tial proton energy spectrum oc E^^-K 

EGRET observations of NGC 253 (Sreekumaret al.T^gTl, 
the nearest spiral galaxy outside the local group (at distance 
d ~ 2.6 Mpc), resulted in an upper limit on the integral flux 
above 100 MeV of 1.0 x 10"^ photons cm"^ s"'. A reanaly- 
sis of these data by Paglione et al. ( 1996) yielded a consistent 
limit of 8 X 10"^ photons cm"^ s~'. More complete EGRET 
data led Blom et al. (1999) to revise this upper limit down to 
3.4 X 10"^ photons cm"^ s"'; for M82 they derive an upper 
limit of 4.4 x 10"* photons cm"^ s Similar upper limits for 
both NGC 253 and M 82 were obtained in a recent reanaly- 
sis of the EGRET data by CiUis et al. ( l2005l . NGC 253 was 
detected in the hard X-ray band using the OSSE instrument 
(Bhattacharya et al. '1994). From the 4.4cr significance detec- 
tion a flux of 3 X 10"" erg cm"^ s"' (50-165 keV) was derived. 

At TeV energies, observations of bofli NGC 253 and M 82 
were made by the HEGRA collaboration, but neither object was 
detected (Getting 2005 1. For NGC 253 a flux upper limit of 
Fy(> 5.2 TeV) < 1.3 x 10"'^ photons cm"^ s"' at 99% con- 
fidence was derived from 32.5 hours of live time. For M 82 
the corresponding limit was Fy(> 2.1 TeV) < 2.7 x 10"'^ 
photons cm"^ s"' from 43.9 hours of live time. The report of 
very high energy (VHE, used here as 0.1 TeV - 10 TeV) y- 
ray emission from NGC 253 by the CANGAROO collabora- 
tion (Itoh et al. I2002> marked the first claim of y-ray emis- 
sion from a starburst galaxy. The differential energy spectrum 
(Itoh et al. 2003b.l can be described by: dF/dE = (2.85 + 
0.71) X 10-12 X {E/lTeY)-^-^^^^'^^ cm^^ s"' TeV"', with a 
corresponding flux above 400 GeV of w 1 .4 x 10"' ' cm"^ s"' , 
or « 15% of the flux of the Crab Nebula above the same thresh- 
old. The emission reported by CANGAROO is extended with 
an rms angular size estimated as 0.3-0.6°, somewhat larger 
than the optical size of the galaxy (28' x 7'). 

Further evidence for particle acceleration in NGC 253, 
at lower photon energies, exists mainly in the form of ra- 
dio data. Radio synchrotron observations show an extended. 
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quasi-spherical halo, rising to at least 9 kpc above the disk 
(Carilli et al. 1992), with a morphology that suggests a gas 
outflow from the disk, carrying with it magnetic fields and 
radiating electrons. Such an outflow should be characteristic 
for starburst galaxies. It has been studied in detail in thermal 
soft X-rays (e.g. Dahlem et al. HWI Strickland et al. EUUni 
I2002al I2002b> and interpreted in terms of "superwinds". The 
observed X-ray emission shows a good correlation with optical 
Ha line and radio continuum observations over the same region. 
From general considerations of particle acceleration, it seems 
likely that accelerated nuclear particles are also present and 
contribute to driving the wind. The CANGAROO detection has 
been interpreted as inverse Compton emission, implying the ex- 
istence of an extended halo of multi-TeV CR electrons around 
NGC 253 (Itoh et al. 2003 ). A cold dark matter (CDM) annihi- 
lation interpretation was also discussed (Enomoto et al. l2003> . 

The detection of NGC 253 by CANGAROO, along with 
the theoretical expectation of very high energy y-rays from 
starburst galaxies, motivated the observation of this ob- 
ject with the High Energy Stereoscopic System (H.E.S.S.). 
H.E.S.S.is an array of four 13 m diameter telescopes, 
each equipped with 960 pixel (5° field-of-view) cameras 
(Hinton 2004). Commissioning of the array in Namibia 
(at 1800 m above sea level) began in 2002 and the ar- 
ray was completed in December 2003. The optical sys- 
tem of the H.E.S.S. telescopes is described in detail by 
Bernlohr et al. ( I2003> and Cornils et al. (20TF3 ), the trigger 
system by Funk et al. (|2QQ4) and the camera electronics by 
Vincent et al. J2003> . Observations of NGC 253, made during 
2003 with H.E.S.S. , are presented here. 

2. H.E.S.S. Observations 

Observations of NGC 253 were made during the construction 
of the H.E.S.S. array. In August and September 2003, 30 hours 
of data were taken with two telescopes in operation. In October 
2003 a further 8 hours of 3-telescope data were taken. After 
run quality selection (primarily to discard data compromised by 
unstable weather conditions) and dead time correction, datasets 
of 24.6 hours (August/September) and 3.4 hours (October) re- 
main. The mean zenith angle of observations was 14°. Only 
events where at least two telescopes provided shower images 
were used in this analysis to enable stereoscopic reconstruc- 
tion. Data were taken in wobble mode, with the target source 
off'set by 0.5° from the tracking position of the telescopes. The 
energy threshold for this dataset is 190 GeV. Two independent 
analysis methods have been applied to search for y-ray emis- 
sion in this dataset. 

2.1. Analysis Method 1 

The standard calibration and analysis used for H.E.S.S. 
is described in detail in (Aharonian et al. 120051 
Aharonian et al. I2004> and is referred to here as Analysis 1. 
Following calibration, tail-cuts image cleaning and image 
moment analysis (Hillas parameter determination) were 
used to characterise each telescope image. Shower direction 
reconsti'uction was made based on the intersection of the major 
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axes of the individual telescope images. The angular resolution 
(defined by the radius containing 68% of the reconstructed 
events from a point-source) is 0.14° for the two telescope 
data set and 0.12° for three telescopes. The primary y-ray 
energy is estimated (with a resolution of ^ 15%) using the total 
image size in each telescope and the reconstructed shower core 
position. The mean reduced scaled width and length of shower 
images are used to select y-ray candidate events. 

Given that the observations were taken in wobble-mode, 
the position 0.5° from the tracking position on the opposite 
side of the field of view (FOV) can be used for a simultane- 
ous background measurement. Figure^shows the distribution 
of squared angular distance (0^) of candidate y-rays from the 
centre of NGC 253 for the combined 2- and 3-telescope data 
sets. The distributions for the on-source and background mea- 
surements are consistent; no evidence for a y-ray signal is seen. 
The vertical dashed lines in this figure indicate the positions of 
the standard H.E.S.S. y-ray selection cut for point-like sources 
(6 < 0.14°) and a much wider cut (at 0.5°) at the approximate 
extent reported by the CANGAROO collaboration. 

To enable 2-dimensional mapping of the FOV a different 
background model is applied. The upper panel of figure |2 
shows the significance of a point-like excess in the FOV of 
NGC 253 using a ring (of 0.5° radius) around each position to 
provide a background estimate (a correction for the radial ac- 
ceptance function of the instrument is made in this case). The 
distribution of significance within the FOV is consistent with 
random fluctuations of the background. 

2.2. Analysis Method 2 

An independent analysis method for the H.E.S.S. data 
(Lemoine-Goumard et al. 2004a Analysis 2) has been used to 
confirm this result. This analysis employs not only different 
reconstruction and background reduction methods but also an 
independent detector calibration (Aharonian et al. 2004 1. The 
reconstruction method is based on a simple 3D-modeling of 
an electromagnetic shower assuming rotational symmetry. For 
each event eight parameters are reconstructed: the direction of 
the shower in the reference frame of the telescope, the core 
position on the ground, the altitude of shower maximum, the 
total number of emitted Cherenkov photons, and the longitudi- 
nal and transverse standard deviations of the Gaussian distribu- 
tions, referred to as "3D-length" and "3D-width", respectively. 
In this method the 3D-width is the most efficient and robust pa- 
rameter to discriminate y-ray showers from the hadronic back- 
ground. 

The background subtraction is performed independently 
in each bin of the FOV using a maximum UkeUhood method 
(Lemoine-Goumard et al. '2004bi based on the knowledge of 
the 3D-width distributions of y-rays and background events, 
respectively. The y-ray width distribution is determined via 
Monte-Carlo simulations and checked on data from known 
point-like sources; the distribution for background events is ob- 
tained from a FOV in the same zenith angle range with no de- 
tected sources. The dependence of these distributions on the 
bin position in the FOV is taken into account. The method 
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Fig. 1. Top: Angular distribution of y-like images relative to 
the centre of NGC 253 ("On") and relative to the background 
control region ("Off"') for Analysis 1. Events are plotted versus 
the squared angular distance to give equal solid angle in each 
bin. Background curves (histograms) are determined relative 
to to points 1 ° away from the source position. Bottom: as for 
upper figure but using Analysis 2 (described in section IT^ . 



yields two independent sky maps, one for y-rays and the other 
for background events. In contrast to Analysis 1, this method 
makes no a priori assumption of azimuthal symmetry in the 
structure of the background. Furthermore, the y-ray contents 
of different bins are statistically independent. This method was 
applied to the same data sets as for the analysis presented 
above; the significance map for y-rays, shown in the lower 
panel of Figure |2 is again compatible with the absence of a 
signal in the whole FOV. 

2.3. Flux Upper Limits 

Figure|3]shows integral flux upper limits calculated assuming a 
spectrum of photon index -3.85 (as reported by CANGAROO) 
and following the method of Feldman & Cousins (il998j . Flux 
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limits are shown for point-like and extended emission (0.5° 
radius) for Analysis 2. Also shown are CANGAROO-II inte- 
gral data points derived from the differential spectrum given 
by Itoh et al. (I2003W. Above 300 GeV we derive 99% confi- 
dence flux upper limits of 1.9 x 10"'^ photons cm"^ s"' for a 
point-like source and 6.3 x 10"'^ photons cm"^ s"' for a source 
of radius 0.5°. The upper limit curves derived from Analysis 1 
are rather similar in shape and have slightly higher values at 
300 GeV: 2.2 x 10"'^ photons cm"^ s"' photons cm"^ s"' 
(point source) and 6.9 x 10"'^ photons cm"^ s"' (0.5°). 

The non-detection of this object cannot be attributed to 
a mis-pointing of the instrument. The systematic error in the 
pointing error of H.E.S.S. during 2003 was < 1' (Gillessen et 
al.|2003 1. Moreover, the correct pointing of the instrument for 
this particular dataset has been independently confirmed using 
the currents induced by stars in the FOV in the photomultiplier- 
tube camera. The limit on the mis-pointing derived from this 
method is < 3'. 



3. Discussion 

Given the strong signal from NGC 253 reported by the 
CANGAROO collaboration, the non-detection of this object 
with H.E.S.S. is surprising. As the source reported by Itoh et al. 
J2003bt is extended in nature, time-variability of the y-ray 
emission cannot be invoked to solve the apparent contradic- 
tion between the H.E.S.S. and the CANGAROO results. In this 
situation it is instructive to review the expectations for y-ray 
emission from this object. 

The speculative possibility of detectable y-ray emission 
from Dark Matter annihilation in the halo of NGC 253 is dis- 
cussed in detail in Enomoto et al. (2003). We concur with the 
view of these authors that such a scenario cannot produce VHE 
y-rays with fluxes detectable with current instruments. The 
standard picture of y-ray production via CR interactions seems 
more plausible given the rapidly growing knowledge about this 
object in other wavelength ranges, principly in terms of super- 
nova explosion rates and the total gas mass in the starburst re- 
gion. The relevant CR interactions are: i) inelastic collisions of 
relativistic nuclei with thermal gas nuclei producing n'* mesons 
which decay into y-rays ii) Inverse Compton (IC) collisions 
of relativistic electrons with ambient low-energy photons iii) 
Bremstrahlung collisions of such electrons with gas nuclei. The 
y-ray flux from these processes is proportional to the density of 
these relativistic particles (CRs). The CR density is determined 
by the strength of their accelerators in the SB region, their prop- 
agation modes, and their energy losses. In particular, propaga- 
tion of CRs via convection in the SB wind is a key process 
that is usually not considered for our own galaxy. In addition, 
the nature of particle diffusion in the SB is likely modified by 
non-linear effects induced by the CRs themselves. 

A self-consistant picture for the y-ray production is dis- 
cussed below. It turns out that the strong localisation of CR 
production and the high gas mass of the SB region lead to a 
rather simple prediction for the VHE y-ray flux, independant 
of the further details of the SB configuration. 



3. 1 . Starburst parameters for NGC 253 

The enhanced CR density and the very high gas and photon 
densities in the starburst (SB) regions are related to each other: 
the high gas densities are necessary for the high star formation 
rate which produces massive stars of high photon luminosity. 
Their subsequent supernova (SN) explosions heat the gas and 
generate CRs. Together the heating of the gas and the high CR 
density drive a mass outflow in the form of a galactic wind. 

NGC 253 has a central SB region which is very bright 
in radio, far infrared, and X-ray frequencies. Following 
Weaver et al. (2002^ we assume it to be a cylinder of ra- 
dius R = 150 pc and height H = 60 pc in the plane 
of the galaxy, with an outflow that is collimated by a mas- 
sive molecular torus into the vertical direction, perpendicu- 
lar to the disk. The outflow velocity Vwind is not very well 
known, with "minimum reasonable values between 300 and 
600 kms (Strickland et al.'2002aJ. However, values for the 
outflowing SN ejecta as high as 3000 kms"' are considered 
possible. For the purposes of our estimates we assume here 
Vwind = 500 kms-'. 

The volume of the SB region is then Vsb ~ 4.2 x 10'' pc^*. 
The SN rate in this small nuclear volume is estimated as 
vsN - 0.03 yr"', about as high as that of our entire Galaxy 
(cf. Engelbracht et al. [^98 1, but values as high as 0.08 yr"' 
have also been estimated (van Buren & Greenhouse 119941 
Paglione et al. 1996i. The rest of the galaxy is thought to 
have a SN rate of the same magnitude, distributed over the 
disk. The total (dominantly molecular) mass M in the nu- 
clear region is rather uncertain: 3 x 10^ < M/Mq < 1.4 x 
10^ (Mauersbergeretal. 1996 Engelbracht et al. 1998 1. We as- 
sume M = 6 X 10^ Mq which corresponds to the preferred value 
of Engelbracht et al. J1998> and gives a density of 580 pro- 
tons cm"^ in the nuclear region. 

3.2. Transport of VHE particles 

The non-thermal properties of such an SB system are quite 
complex, as the very extended radio synchrotron emission 
shows. This radio emission is intimately connected to the large 
expanding halo in thermal X-rays and optical H„-line photons. 
Nevertheless, the TeV y-ray emission can be discussed in a 
fairly straightforward manner. For this purpose, let us assume 
that particle diffusion can be calculated from the theory of the 
wind from our Galaxy (Ptuskin et al. 1997 1, scaled to the SB 
parameters of NGC 253. We point out that in the theory of 
Ptuskin et al. the diffusion coefficient in the Galactic wind is 
assumed to be due to Alfven waves, excited by the particles 
themselves, as they stream away from their sources in the disk. 
The wave amplitudes are limited by nonlinear wave damping 
on the thermal wind plasma. 

To estimate the spatial extent and magnitude of the y- 
ray emission as a result of the SB, we first consider the 
time scales for energetic particle transport and their colli- 
sional/radiative energy losses. For protons the convection time 
tconv - H/Vv/ind ~ 1.2 X lO*" yr, whereas the loss time in the 
nuclear region due to inelastic collisions with nuclei of the gas 
amounts to fioss = 8.6 x lO'' yr. 
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The diffusion vertically to the disk is given by the diffu- 
sion coefficient *■ = (fi/BgaOCe/egaO^'fgaiC^^), where above a 
few GeV k„^i(E) ^9x 10^"^ (fi/lTeV)'' cm^s"' is consistent 
with the mean matter thickness traversed by the CRs in our 
Galaxy (and its rigidity dependance), before they are ultimately 
convected away into extragalactic space (Ptuskin et al. 119971 
In this scaling relation e oc vsN/fC/? source surface area) de- 
notes the spatial energy flux density of CRs released into the 

where is the 

CR source power density, as a function of particle momentum 
p, per unit surface area of the source region. As we shall see, in 
this case fioss is smaller than both fconv and fdiif (which is defined 
below). This results in an upper limit for the y-ray emission 
from NGC 253 for a given energy input into CRs. 

To determine e for NGC 253 we note that the cylindrical SB 
region has a radius of 150 pc, according to the above model, 
whereas the corresponding region of our Galaxy, the disk, has 
a radius of 15 kpc. Therefore, assuming vsn = 0.03 yr ' in 
the SB as in the Milky Way, yields e/fgai - 10"*. For the 
magnetic field B in the SB region we chose B < 270 juG, 
where 270 /iG corresponds to the equipartition value derived 
by Weaver et al. t.2002 ) which we consider the maximum pos- 
sible field in the SB region. With a Galactic field strength 
figai = 3 juG this massive CR production reduces the diffusion 
coefficient to k < 9 x lO^^^/fgai. The corresponding diffusion 
time fdiff = H^/k > 1.3 x 10^ yr for a 1 TeV particle also 
exceeds fioss- This means that CR protons with energies up to 
about 1 TeV behave calorimetrically: i.e. they lose all their en- 
ergy by inelastic collisions before being able to leave the SB re- 
gion. However, the increase of k with particle energy makes the 
particles with energies above several TeV diffusion dominated, 
which implies that they escape the SB region before interact- 
ing. The turnover point is at y-ray energies of a few hundred 
GeV. 

For TeV electrons the loss time is estimated as follows: the 
enormous far infrared (FIR) SB-luminosity Lsb = 1.1 x IO^'L© 
deduced by Engelbracht et al. ( 1998 ) implies a radiation energy 
density U^^, ~ 500 eVcm^^*. This leads to an inverse Compton 



(IC) loss time f 



ic 



180(£'/4TeV)-' yr for the electrons of 



energy E ^ 4 TeV that produce 1 TeV y-rays. The IC scattering 
is here considered in the Thomson limit which applies for FIR 
target photons. Therefore, comparing with the proton emission, 
the SB region acts even more calorimeterically for the emission 
of TeV y-rays via IC scattering. 

3.3. VHE y-ray emission 

The hadronic y-ray emission in this scenario simply corre- 
sponds to the total amount of CR energy created. The y-ray 
energy spectrum corresponds to that of the SNR sources, as- 



' In our galaxy V„ind oc ^ for >s <15 kpc where s is the height 
above the disk, The distance s^, within which the energetic particles 
can diffusively return to the disk to produce secondaries is given by 
K^.^](E)/s^ = V„in(j(j,t)/i,t (the diffusion-convection boundary). Thus 
S-i, oc kI^^{E) oc E^-^^ for our galaxy. The probability for a particle to 
cross the disk, and thus the mean matter thickness traversed by parti- 
cles of energy E is inversely proportional to s^, . 



sumed here to be dNcR/dE oc E^^ \ In this picture diffusion is 
irrelevant below a few TeV, and the total CR energy Eqr in the 
SB region amounts to 



£cR* 1.5x1 0^^ erg 



0£sN 



teif 



0.03yr-i lO^Oerg 5 x lO^yr 



(1) 



where f , 



= tu 



+ L ' „ + t,L, and < 1 is the CR en- 



efF loss 'conv ^ 'diff' 

ergy production efficiency for an average SN event with a me- 
chanical energy release Esn- The CR energy density E^^^ - 
/iCR/Volume amounts to E^^^ ^ 650 eVcm"-'. This has to be 
compared to the Galactic value Ef^ ~ 1 eVcm"-'. Assuming 
the radiating TeV particles to freely penetrate the gas in and 
around their sources - a premise which would become more 
and more questionable at lower energies on account of the 
shorter scattering mean free path A - 3k{E)/c - the energy 
spectrum of the hadronic y-ray flux can be calculated from the 
following relation (cf. Volk et al. 1996 1: 



Fy(> E) = 1.7 X 10"'^ photons cm"^ s" 



UTeV/ \leVcm-3/\lMpc/ 



M 



\ lO'^Mp 



(2) 



Substituting our preferred values for d, M, and E^, as cal- 
culated above, we obtain Fy(> E) = 9.8 x 10"'^ x 6 x 
(/i/lTeV) photons cm ' s ' up to a few hundred GeV, and 
significantly falling off at higher y-ray energies. The factor 6 
denotes the uncertainty of the assumed parameters, mainly vsn, 
M, and 0. Assuming that » 0.1 and the average /isN - 
10^' erg, we obtain 0.5 < 6 < 5. 

As long as M > 6 X IO^Mq we have feff ~ floss °^ Vsb/M. 
Since E^ oz fioss/VsB, Fy (where Fj.cc E^ ■ M) then becomes 
independent of M, Vsb, fconv and fdiit, and is only proportional 
to the energy input per SN explosion ©/ssn and the SN rate vsn, 
apart from the geometrical factor c/"^. This is the consequence 
of the calorimetric behaviour of the system. 

The CR electron energy produced in the SNR sources of 
the SB can be assumed to amount to no more than a fraction 
~ 10"^ of the energy generated in nuclear CR particles, as is 
the case in our galaxy. The contribution of IC scattering to the 
y-ray emission in the VHE region can therefore be safely disre- 
garded. The above hadronic flux value, with ^ = 1 , is similar to - 
in fact somewhat larger than - the H.E.S.S. upper limit for the 
VHE y-ray flux. At the most relevant energy of about 200 GeV 
the predicted flux is about a factor of 3 larger than the H.E.S.S. 
limit for an assumed point source. This demonstrates the phys- 
ical relevance of the H.E.S.S. upper limit measurement. 

Our discussion rests on the assumption that the propaga- 
tion properties of energetic particles in our own Galaxy and its 
halo can be scaled to the SB parameters of NGC 253. If this 
assumption does not hold, then some fraction of the acceler- 
ated particles will escape the galaxy before radiating, and the 
predicted y-ray flux from NGC 253 will decrease relative to the 
estimate given here. 

The emission from the rest of the galaxy may be spatially 
extended. However, assuming a situation roughly equal to that 
of the Milky Way, where the amount of matter traversed by 
a TeV particle corresponds to less than 1.5% of the hadronic 
interaction length, the majority of the protons will escape and 
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will not transfer their energy to y-rays. Therefore their total 
y-ray emission is substantially smaller than that from the SB 
region. The IC y-ray flux from the associated electrons is still 
negligible compared to the CANGAROO flux. 

This scenario leads to the prediction that VHE y-ray emis- 
sion from the SB in NGC 253 will be highly localised, with an 
extent of less than 30 arcseconds. Such emission would appear 
point-like for current VHE y-ray instruments. As discussed ear- 
lier, extended emission from the galaxy as a whole is expected 
to be rather faint in comparison to the flux reported by the 
CANGAROO collaboration. In addition the expected spectrum 
should be rather hard as discussed earlier, again in contrast to 
the reported CANGAROO spectrum. 



4. Conclusions 

The recently completed H.E.S.S. instrument is currently the 
most sensitive detector of y-rays at VHE energies. The non- 
detection of NGC 253 by H.E.S.S. is surprising given the rather 
high flux previously claimed by Itoh et al. ( I2002> . From theo- 
retical considerations of y-ray production in starburst galaxies 
it appears that the H.E.S.S. limit is close to the expected flux 
from the NGC 253 starburst, already excluding extreme param- 
eter values (6 > 1). However, a diffusive transport away from 
the starburst region that is faster than assumed here cannot be 
excluded experimentally. The possible y-ray flux would then 
be correspondingly reduced. 

It also seems likely that VHE y-ray emission from 
NGC 253 should be rather localised (indeed point-like for cur- 
rent detectors), given the characteristics of the starburst mea- 
sured at other wavelengths. A hard y-ray photon index (~ 2) 
would also be expected. 

Despite this non-detection, NGC 253 remains one of the 
most promising starburst galaxies for detection in very high 
energy y-rays. The likely calorimetric character of this system 
makes deeper observations an interesting proposition. 
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Fig. 3. Upper limits from H.E.S.S. on the integral flux of y- 
rays from NGC 253 (99 % confidence level) from Analysis 
2. Curves for point like emission (solid lines) and for a 
source of 0.5° radius (dashed lines) are shown. CANGAROO 
integral data points (derived from the differential spectrum 
given by Itoh et al. (,2003b^) and an upper limit from 
HEGRA ('Gotting [2005> are shown for comparison. 



Fig. 2. Top: y-ray point source significance map (grey-scale) 
derived from two telescope H.E.S.S. data taken in August 
and September 2004 and using Analysis 1. The white con- 
tours show the optical emission from the galaxy (in a lin- 
ear scale). These data are taken from the STScI Digitized 
Sky Survey 2004 Black contours are confidence levels (at 
40%, 65% and 80%) for TeV y-ray emission reported by 
Itoh et al. (i2002i) . The dashed white line shows the angular cut 
used to derive the extended source flux limits presented here. 
In this approach bins are correlated within the radius of the 
cut (0.14°). Bottom: Significance map derived from the same 
dataset using Analysis 2. In contrast to the Top panel the sta- 
tistical significance in each bin is independent. The bin size is 
matched to the angular resolution of the instrument for these 
data. 



